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Abstract 

Cross sections and analyzing powers for proton elastic scattering from 116, 118,120, 122, 124gj^ 295 
MeV have been measured for a momentum transfer of up to about 3.5 fm^^ to deduce systematic 
changes of the neutron density distribution. We tuned the relativistic Love-Franey interaction to 
explain the proton elastic scattering of a nucleus whose density distribution is well known. Then, 
we applied this interaction to deduce the neutron density distributions of tin isotopes. The result of 
our analysis shows the clear systematic behavior of a gradual increase in the neutron skin thickness 
of tin isotopes with mass number. 
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I. INTRODUCTION 



Charge distributions in stable nuclei have been reliably measured by electron elastic 
scattering and muonic x-ray data [1]. These charge- sensitive experiments have provided 
precise information on charge distributions. On the other hand, it is much more difficult 
to deduce neutron density distributions, since electromagnetic interaction provides little 
information on neutron density distributions. The proton and neutron density distributions 
have a similar shape in stable nuclei. However in recent research it has been reported that 
in some unstable nuclei the differences between proton and neutron shapes are greater than 
those in stable nuclei m. It has also been indicated that the thickness of the neutron skin 

fi n 

is closely related to the symmetry term of the equation of state (EOS) [3|, Thus, the 
determination of neutron density distributions has become increasingly important. 

There have been many experiments attempting to extract neutron and matter density 
distributions in the nuclear interior using hadronic probes [5|]. Pion and alpha elastic sea ;- 
tering have been performed in the study of neutron and matter density distributions 6|, [Zl, l8| . 
Compared with other hadronic probes, the elastic scattering of protons at intermediate en- 
ergies is suitable for extracting information on the nuclear surface and interior, because at 
intermediate energies, proton elastic scattering has a simple reaction mechanism. To deduce 
nuclear densities using protons, the incident energy has to be sufficiently high to describe 
the scattering by the simple mechanism. At energies above 100 MeV, we can explain pro- 
ton elastic scattering microscopically, because the imaginary part of the optical potential 
describes a quasi-free process mainly without the need for a renormalization factor. So far, 
energies above 500 MeV have been applied for proton elastic scattering to study neutron 
density distributions j^, Q. However, this energy is sufficiently high to produce mesons, 
and information on the nuclear interior is easily masked by the meson-productions. Further- 
more, the total cross section of nucleon-nucleon scattering shows a minimum at the incident 
energy of 300 MeV. We thus adopt 300 MeV protons in this work as probes for information 
on the nuclear interior. 

In our previous papers [llj, we tuned the effective relativistic Love-Franey interaction 
by introducing 'so-called medium effects' for the scattering from a nucleus whose density 
distribution is well known. We used elastic scattering from ^^Ni to tune the interaction, since 
^^Ni is the heaviest stable nucleus with N ^ Z and the density distribution of the neutrons 



2 



in ^^Ni can be assumed to be the same as that for the protons. To explain the results of 
our experiments we found that we had to modify the scattering amplitudes of the nucleon- 
nucleon interactions inside the nucleus as follows. We phenomenologically changed the 
masses and coupling constants of exchanged mesons depending on the nuclear density. Thus, 
we could explain the scattering sufficiently well to deduce the matter density distribution 
precisely. 

For our first systematic search for neutron density distributions, we selected tin isotopes. 
Tin has many stable isotopes (^^^Sn-^^'^Sn). Also, unstable tin isotopes have a long isotopic 
chain including two double-magic nuclei (^°°Sn [N=50], ^^^Sn [N=82]). Moreover, its proton 
number is a magic number (Z=50). Thus, tin isotopes are suitable for the study of systematic 
changes in neutron density distributions. The main purposes of this work are to attempt 
to deduce information on neutron density distributions, and to systematically study the 
neutron skin thickness of tin isotopes. The experimental setup and procedure are described 
in Sec. II. Details of the analysis used to deduce neutron density distributions are given in 
Sec. III. The deduced radii of tin isotopes are discussed in Sec. IV. A summary is given in 
Sec. V. 



II. EXPERIMENT 



The measurements were performed at the Research Center for Nuclear Physics (RCNP), 
Osaka University. Polarized protons from a high-intensity polarized ion source were injected 
into an AVF cyclotron, transported to a six-sector ring cyclotron and accelerated up to 295 
MeV. The polarization axis was in the vertical direction. The spin direction and magni- 
tude of the beam polarization were measured continuously using sampling-type beam-line 
polarimeters (BLPs) 1^ placed between the ring cyclotron and a scattering chamber. Each 
polarimeter utilized left-right asymmetries in p-H scattering from (CH2)n foil to determine 
the vertical transverse component Py of the beam polarization. The typical beam polariza- 
tion was 65%. Then the beam was transported to a target center in the scattering chamber. 
The typical beam spot size on the target during measurements was 1 mm in diameter. Fi- 
nally, the beam was stopped by an internal Faraday cup inside the scattering chamber in the 
case of forward-angle measurements. In the measurements at backward scattering angles, 
the beam was transported to another Faraday cup located inside the shielding wall of the 
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experimental room about 25 m downstream of the scattering chamber. The integrated beam 
current was monitored using a current digitizer (Model lOOOC) made by BIG (Brookhaven 
Instruments Corporation). Additionally, the beam current was monitored independently 
using p-p cross sections at the BLPs during the backward-angle measurements. Five tin 
isotope targets (H^, 118,120, 122, i24gj^-j ^^le form of self-supporting metal foils were used for 
this experiment. Two different thicknesses were used for each target. Thin targets were 
used for the forward-angle measurements to reduce the dead time of the data acquisition 
system, and thick targets were for the backward-angle measurements to increase the yields. 
The enrichment and thicknesses of each target are shown in Table [H 

TABLE I: Target enrichment and thicknesses of tin isotopes. 



Nucleus Enrichment Thin Thick 

ii^Sn 95.5% 10.0mg/cm2 lOO.mg/cm^ 

i^^Sn 95.8% lO.Omg/cm^ lOO.mg/cm^ 

i20Sn 98.4% 5.12mg/cm2 39.9mg/cm2 

i22Sn 93.6% 10.5mg/cm2 85.4mg/cm2 

i24Sn 95.5% 5.00mg/cm2 62.7mg/cm2 



The main contaminants of the targets originated from other tin isotopes. The present 
energy resolution could not separate the elastic scattering of other tin isotopes. Thus, we an- 
alyzed the targets including the contamination from other isotopes at all momentum transfer 
regions. We estimated the error of this analysis to be less than 1% for all cross sections and 
analyzing powers. We used an automatic target changer system in this experiment to reduce 
the systematic errors of relative cross sections between isotopes. We formed a stack of three 
targets, which were moved vertically every 2 min to avoid errors due to the drift of the beam 
direction and that of its position on the targets. 

We used a high-resolution (p/Ap ~37,000) magnetic spectrometer, 'Grand Raiden' (GR), 
which had a Q1-SX-Q2-D1-MX-D2 configuration and focal-plane counters for momentum 
analysis p^]. The focal-plane counters consisted of two vertical-drift-type multi wire drift 
chambers (VDGs) and two plastic scintillating counters. The momentum of protons scattered 
from the target was analyzed using the GR. The trajectory of the scattered protons from the 
target was reconstructed from the position measurement carried out using two sets (XI, Ul 
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FIG. 1: Sample focal plane spectra corresponding to the excitation energy of Sn, taken at a 
laboratory scattering angle of 35.5° 

and X2,U2) of VDCs, which were placed near the focal plane of the spectrometer jl^ . 
Each set had an effective area of 120 cm [width] x 10 cm [height]. Two 1 cm-thick plastic 
scintillating counters (PSl and PS2) placed behind the VDCs were used as triggers and for 
particle identification. 

The energy resolution of the beam was 200 keV in FWHM, which was determined by the 
energy width of the beam itself and was sufficient to separate elastic peaks from inelastic 
peaks. Figure [1] shows a sample focal-plane spectra of ^^°Sn. The differential cross sections 
and analyzing powers were measured for analysis of up to 50°, corresponding to a momentum 
transfer of 3.5 fm~^ in the center-of-mass system. We included 3% errors for both cross 
sections and analyzing powers as errors resulting from experimental conditions including 
target thickness uncertainty. 

III. THEORETICAL ANALYSIS 

We analyzed our data using the formula for the relativistic impulse approximation (RIA) 
using the relativistic Love-Franey interaction [15]. The invariant amplitude has been deter- 
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mined from nucleon-nucleon phase shifts and is expressed as 

F{q) = F' + F^5,)7[;) + F^'lMi) 

+^Xo?(i)M^ + ^^(o)7ro)7fi)7(i).. (1) 

Each amphtude is shown as the sum of real and imaginary amphtudes using the masses, 
couphng constants, and cutoff parameters of exchanged mesons. Figure [2] shows the exper- 
imental results for five tin isotopes and two different theoretical results. The solid lines in 
Fig. [2] are the results of RIA calculations with the default parameters proposed by Murdock 
and Horowitz 



ISjand the relativistic mean fielc 



the calculations using the global potential |l7 



RMF) densities jl6|]. The dashed lines are 



181] . While both calculations are in good 



agreement with the experimental data of obtained from the analyzing powers, the differ- 
ential cross sections are overestimated, particularly in the large-momentum-transfer region. 
Murdock and Horowitz used different masses and coupling constants for real and imaginary 
scattering amplitudes, depending on the interaction. In our previous work, we attempted 
to explain the experimental data by phenomenologically changing the masses and coupling 
constants of exchanged mesons (a, uj) in the real and the imaginary amplitudes depending 
on the nuclear density. The formula for density dependence is as follows; 

2 -2 9j 9j 

9i,9 



mj,mj — ^ mj[l + bj p{r)/po],mj[l + bj p{r)/po], (2) 

where, mj, rhj, gj, and, gj indicate the masses, and coupling constants of mesons for real and 
imaginary amplitudes, respectively, where j refers to the a, uj mesons. The normal density 



Po is 0.1934 fm ^ Uj, ll9l]. These changes in the masses and coupling constants are called 



medium effects, and may be an effect of the presentations of the partial restoration of chiral 
symmetry, Pauh blocking, and multistep processes. The tuned effective interaction is applied 
to existing ^°^Pb data obtained at TRIUMF [2^ with a density distribution calculated from 
the RMF in Fig. [31 It was found that the tuning of the effective interaction was meaningfully 
improved compared with the original unmodified interaction. 

For the RIA calculation, we need four density distributions; the vector and scalar density 
distributions of protons and neutrons. The proton density distribution and the relation 
between the scalar and vector density can be obtained from the charge distributions, the 
nucleon electric form factors, and the RMF calculation. Thus, we can determine the neutron 




FIG. 2: Differential cross sections and analyzing powers for proton elastic scattering from tin 



isotopes. The solid lines show the results of RIA calculations with the default parameters proposed 



by Murdock and Horowitz 
using the global potential 
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and the RMF densities [3], while dashed lines show calculations 
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FIG. 3: Experimental data for ^"^Pb measured at TRIUMF [20|lat 300 MeV. Dashed lines show the 



result of the original RIA calculations using RMF densities 



16l |. Solid lines show the calculation 



based on the modified effective interaction using RMF densities. 



density distribution by comparing experimental data with a calculation using the tuned 
effective interaction. 

In this paper, we attempt to extract neutron density distributions for tin isotopes, con- 
sidering the various ambiguities caused by the modification parameters used in the RIA 
calculation, the proton form factors, and an assumption based on scalar densities. 



A. Proton density distributions of tin isotopes 



For the RIA calculations we used point proton density distributions derived from charge 
distributions observed in electron-scattering experiments [l]. We used the sum-of- Gaussian- 
type (SOG) charge distributions of ^^^'^^^Sn, which were obtained from the electron scat- 
tering over a 3.5-fm~^-wide momentum transfer region The charge distributions are 
described as 

Z ^ Q, 



Pch{r) 
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FIG. 4: Existing experimental results for charge radius of tin isotopes. Squares and triangles 
represent results from electron elastic scattering using different shapes of charge distributions l|. 
In the case of SOG data, from existing obtained -"^-^^'^^^Sn data, 118.120,122 g^^ data are deduced in 
this work. Circles show results from muonic x-ray data [23]. Diamonds show combined results 



from both electron scattering and muonic x-ray |22l | 

However, reported data on charge distributions obtained using model-independent distri- 
butions are limited. Figure H] shows existing experimental results for the charge radius of 
tin isotopes. Since the radius of tin isotopes increases smoothly with mass number and the 
number of protons is constant, we can expect a smooth change in the charge distributions 
as the neutron number increases. Thus we obtained the charge distributions of ii8.i20,i22gj^ 
by interpolation using the arithmetic mean of the previous derived charge distributions of 

116,124Sn; 



p::'» = ""-^ "^^'"-^ (4) 



n^^''Sn(^'. _ Pch > + Pch fK\ 
Pch V) — 2 ^ 



Pch {' ) — I • y^j 

The interpolated radii of model-independent-type charge distributions for ii8,i20,i22gj^ 
consistent with other experimental data obtained from the electron scattering |l| and muonic 



x-rays [23|, as shown in Fig. |H The differences between interpolated radii and other results 
from the combined analyses of elastic electron scattering and muonic x-ray data elastic 
electron scattering and muonic x-ray data are small 



221 ] ■ We estimate that the errors of our 
deduced charge radii are 0.003 fm, which are the same order as those for ^^^'^^^Sn 

To extract the point proton density distribution we need to unfold the charge distribution 
with the finite size of the protons and to consider the contribution from neutrons. Therefore, 
we can write a charge distribution pch using the charge distributions of a proton p^^°*°" and 
neutron as follows; 



pch{r) = j pp{?)f^j,^°*'""{f-?)dr' 

+ [ pr,{P')p2'''""'{r-r~^')dr^'. (7) 



Therefore, we calculate the mean square radius of the charge distribution as follows; 
where 



(rlh) = j r^pch{r)df/ j Pch{r)dr 
(rl) = J r^Pp{r)df/ j Pp{r)df 
{{vT^y) = / r'pT"''{^dr/ j pT'^if^df 



('^ch)^^^ ' ('^p)^''^ ' (('^c™*°")^)^''^ ' ^^"^ ((^ch"*™")^)^'^^ denoted as the root-mean-square 
(RMS) radii of pch, Pp, and p^ewtron^ respectively. We used the simple parameteriza- 

2J], where the RMS radius of the proton itself is 0.863(4) fm 



tion of nucleon form factors 



and the mean square radius of the neutron itself is -0.112(3) fm^. This proton radius is con 



sistent with that obtained from the measurement of the hydrogen IS Lamb shift [25|]. The 
contribution to the RMS radius of the point proton density distribution from the neutron 
form factor is in the order of 0.02 fm. Thus, the contribution from the neutron form factor 
is small but not negligible. 

Equation ([7]) shows that we need the neutron density distribution to precisely deduce 
density distribution. Since the errors of proton radii are mainly determined by the charge 
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radii in the equation, the errors of the RMS radii of the point proton of tin isotopes are 
estimated to be 0.003 fm. 

In this report we tune newly a medium-effect parameter set using the ^^Ni data of 
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26| . because we adopted a new treatment for the point proton density distri- 



Refs 
bution. 



B. ^^Ni and the effective interaction 

We must cahbrate the effective nucleon-nucleon interaction, particularly in the nuclear 
interior, by the scattering from the nucleus. In this subsection we discuss the analysis of 
the elastic scattering of ^'^Ni at the energy of 295 MeV, which allows to adjust the effective 
interaction further. 

In the previous work, eight parameters were searched for independently including the 



imaginary parts of two exchanged mesons HI]. It is difficult to obtain a unique medium- 



effect parameter set due to its many degrees of freedom. This difficulty causes ambiguities in 
the neutron density distributions. Thus, we attempt to express the elastic scattering using 
the medium-effect parameters that have less freedom. We need the information of both 
proton and neutron density distributions to calibrate the effective interaction based on RIA 
calculations. Several calculations have been performed to obtain the proton and neutron 
density distributions of ^^Ni using relativistic or nonrelativistic mean-field calculations. The 
results are dependent on the interactions used in the calculations but the differences between 



the RMS radii of point protons and neutrons, the neutron skin, are generally small 
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271, 



28|, 



32], and are given as 

Ar.,^{riy^'-{riy^\ (9) 
The used of SOOMeV proton elastic scattering based on KMT nonrelativistic impulse ap- 

n 

proximation gave a thickness of ArVp=0.01(5) fm [9|. X-ray measurement of an antiprotonic 



atom gave a value of -0.9(9) fm [33|]. Thus, in the case of ^^Ni, we can assume that the neu- 
tron density distribution has the same shape as the proton density distribution. Therefore, 
the neutron density distribution can be described as 

p„(r) = (iV/Z)p,(r). (10) 



11 



1 



Proton (RMF) 
Neutron (RMF) 
Trial Ratio 

I .... I .... I ... . 

0.01 0.02 0.03 0.04 O.OS 0.06 0.07 0.08 0.09 0.1 

Proton Vector density(f m"') 

FIG. 5: Ratios of scalar density to vector density for proton and neutron calculated by the RMF 
calculation for ^^Ni. 

Thus, the substitution of Eq. ffTOl) into Eq. ([7]) gives Eq. ffTTj) the point proton density 
distribution of ^^Ni, 

Pchir) = j Pp(0[Pcr°"(^-^) 

+ (iV/Z)p;!^^"*"°"(f-r')]rfr'. (11) 

According to the RMF calculation, the ratio of scalar-to-vector densities has an almost 
constant value of 0.96, as shown in Fig. [5|, corresponding to the ratio of the integrated scalar 
density to the integrated vector density. Using this constant value, we reproduced the scalar 
density distributions from the RMF calculation. Even if we had taken different constant 
value from 0.95 to 0.97, the differences in the cross sections and analyzing powers would have 
been less than 2-3%, which can be compensated for by the ambiguity of the medium-effect 
parameters. 

Since in the case of eight parameters, we have too many degrees of freedom for fitting 
parameters and the correlations among these parameters are large, in this work we have 
used the same modification parameters for both real and imaginary parts of the scattering 
amplitude {a.j=a.j ,bj=bj). Figure [6] shows the correlation between modification parameters 
of ao- and a^. The correlation is valley like and very strong; the correlations of the other five 
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combinations of modification parameters also show similar behavior. Thus, we determined 
the modification parameters uniquely. 

The medium-effect parameters were determined by fitting the experimental data so as to 
minimize the chi-square (x^) value, is defined as 

id — [(^^exp. ~ ^theo)/ ^^exp] (12) 

where ^exp., Axexp. and yitheo. are the experimental data, the errors in the data, and the 
calculation results, respectively. Figure [7] shows the experimental data at 295 MeV and the 
fits with the RIA calculation using the modified effective interaction with the unfolded proton 
and neutron densities of ^^Ni obtained using Eq. ffTOj) . Our calculations are in good agreement 
with all the data of the cross sections, the analyzing powers, and the spin rotation parameters 
for up to 3.5 fm~^. Xmm; ^1^6 minimum of for ^^Ni data using four free parameters, has 
almost the same order as that in the previous work using eight free parameters. We estimated 
the errors of the modification parameters from the statistics and the experimental conditions 
as follows. 

< xLn + P (13) 



Here, P (in this work P=4) denotes a number of fitting parameters 3J]. Table HTl shows a 
summary of the fitting results. Figure [6] shows the correlation between a^ and a.i^. The 
circle shows the best-fit parameter, and horizontal and vertical bars show the error bars of 
ao- and a^^, respectively. The hatched area shows the region that satisfies Eq. (|T3l) . Since 
each pair of parameters has a strong correlation, we cannot independently use the entire 
region defined by the errors in Table Ull meaning that the usable area defined by each pair 
of the parameter set is very narrow. 

TABLE II: Medium-effect parameters at 295 MeV. The error estimation is discussed in Sec. II. B 
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FIG. 6: Contour plot of correlation of ao- with a^u for ^^Ni. The hatched area represents the 
region obtained from Eq. (jl3p . The open circle and bars represent the best-fit parameter and the 
errors shown in Table HH respectively. The magnitude of is arbitrary. 

C. Tin isotopes and neutron density distributions 

We have tried to deduce the neutron density distributions of tin isotopes using the newly 
modified effective interaction. The scalar nucleon density distributions of tin isotopes were 
assumed to be 0.96 times the vector nucleon density distributions as in the case of ^^Ni, 
because the ratios of the scalar to vector densities for tin isotopes are almost the same [l^. 
The initial proton density distributions of tin isotopes are obtained by Eq. ([7]) using the 
neutron density distribution given by Eq. ( JTOl) . We used a SOG for the neutron density 
distribution, which is expressed as 



( \ = N sr 



X 



(14) 



Here, N is the number of neutrons. Equation (HM is almost the same as Eq. (I3l). Whil( 
parameters such as 7 and Ri were fixed using results of the charge distribution of ^^^'^^''Sn [l| 
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FIG. 7: Experimental data for ^^Ni at 295 MeV and results of several calculations. Dashed lines 
show results of the original RIA calculations using RMF densities. Dotted lines represent similar 
calculations but using densities deduced from electron-scattering data. Solid lines show calculations 
based on the modified effective interaction using the deduced densities. 



Qj were determined independently under the following normalization condition; 



/12 
p„(r)rfr = iV^^g, = l. 
i=i 



(15) 



For fitting, we used a range of momentum transfer from 0.7 fm^^ to 3.5 fm~^, since the 
data for the forward angle (^c.m. < 8°) were difficult to measure experimentally. Our range 
covers the same momentum transfer region of electron scattering as that in Ref. For 
the medium-effect parameters, we used the same parameters as those determined by the 
scattering from the ^^Ni target. 

In Fig. [HI we show experimental data together with the fitting results and calculations 
using the initial neutron density distributions, which are defined as 

Pn(r) = p,(r) + (pr^(r) - pf ^^(r)), (16) 

where p^f^fif) denote the RMF neutron (proton) density distributions {l^. We determine 
the proton and neutron density distributions by iterating Eq. ([7]) several times until self- 
consistency is achieved. The process of iteration hardly affects the proton densities, because 
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FIG. 8: Differential cross sections and analyzing powers for proton elastic scattering from tin 
isotopes. The figure includes data from Fig. [2l The dashed lines show calculations based on 
the modified effective interaction using trial density distributions. The solid lines are best-fit 
calculations based on the modified effective interaction. 
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the proton density distributions are not sensitive to the fine structure of the neutron density 
distributions. 

The calculations using the initial neutron density and the medium- modified effective 
interaction reproduce the absolute value of the cross sections and analyzing powers. The 
angular distributions of both the differential cross sections and the analyzing powers using 
the initial neutron density are shifted to forward angle in angle because the RMF calculation 
has a tendency of a larger neutron radius. 



D. Uncertainties of neutron density distributions 

We estimated the uncertainties of the neutron density distributions of the tin isotopes. 
There are generally two types of uncertainty. One originates from the experiment, and the 
other originates from the model of calculation. In our case, the reduced chi-square "x^jv 
is larger than 1, and is typically 7 for Ii6,ii8,i22,i24gj^_ ^.j^g ^^^^^ ^^°Sn, not only the 
cross section and analyzing powers but also the spin rotation parameters [26| are used as 
experimental data in the analysis. Thus, the errors for ^^°Sn resulting from the experimental 
conditions, degrees of freedom v and x^/ ^ are different from those for other isotopes, x^/ ^ is 
about 10 for ^^°Sn. If the model were perfect, it would realize j v=\. Thus, the reason that 
'X^ jv >1 is attributed to the inadequacy of our model. To compensate for this inadequacy of 
the model, we increased the errors artificially by multiplying all the experimental errors by 
a constant factor to realize 'y^ lv=\. This means that we must redefine Eq. (fT3!) as follows; 

< xLn + ^ X xLn/'^- (17) 

The uncertainties of the neutron density distributions are estimated by the following proce- 
dure. 

First, we obtain a new medium-effect parameter set using ^^Ni data by the Monte Carlo 
method. We adopt the parameter sets that satisfy Eq. ( fT3|) as 'good' medium-effect param- 
eter sets. 

Second, we estimate the uncertainties of neutron density distributions using each set 
of good medium-effect parameters. The neutron density distributions are calculated by 
the Monte Carlo method under the normalization condition of Eq. (fT5l) . We compare the 
experimental data with the result of the calculation using each good set of medium-effect 
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FIG. 9: Point nucleoli density distributions for tin isotopes. Solid lines show point proton density 
distributions. Dashed lines show best-fit neutron density distributions. Hatched areas represent 
the error envelopes encompassing all the trial SOG neutron density distributions deduced by the 
modified effective interaction. Not only statistical and experimental errors but also systematic 
errors from the model dependence contribute to these regions. 

parameters and the trial neutron density distributions. We also adopt all neutron density 
distributions that satisfy Eq. (fT7|) as good neutron density distributions. 

We repeated this procedure until the uncertainty of each good neutron density distribu- 
tion converged. The uncertainties of the neutron density distributions in this analysis are 
defined as the outskirts of the density distributions at each radial point. The uncertainties of 
the RMS radii of the neutron density distributions were also calculated using all the neutron 
density distributions that satisfy Eq. (ITTI) . 

Figure [9] shows the deduced neutron density distributions and point proton density distri- 
butions of the tin isotopes. Hatched areas in the figure represent outskirt regions encompass- 
ing all the neutron density distributions for tin isotopes allowed by Eq. ( fT7|) . corresponding 
to the uncertainties of the neutron density distributions in this analysis. Table IIIII shows a 
summary of the RMS radii of the proton and neutron density distributions and the neutron 
skin thicknesses of the tin isotopes. 
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TABLE III: Reduced RMS radii and thicknesses of tin isotopes. 



Target (p^)^/' (fm) (p.^)^/' (fm) Ar^p (fm) 

ii'^Sn 4.562 ± 0.003 4.672 ± 0.018 0.110 ± 0.018 

ii^Sn 4.575 ± 0.003 4.720 ± 0.016 0.145 ± 0.016 

i20Sn 4.589 ± 0.003 4.736 ± 0.033 0.147 ± 0.033 

i22Sn 4.602 ± 0.003 4.748 ± 0.016 0.146 ± 0.016 

i24Sn 4.615 ± 0.003 4.800 ± 0.017 0.185 ± 0.017 



In the high-momentum-transfer region, the multistep process may affect the angular 
distribution of elastic scattering [35|. To estimate the effect of multistep processes, we 
performed a coupled-channel calculation by using the coupled-channel code ECIS95 written 
by Raynal jsol . Figure [10] shows the calculated cross sections and analyzing powers with 
and without performing a coupled-channel calculation using the global potential 17, isl . 
The contribution from the coupled-channel calculation for ^^Ni in the momentum transfer 
range from 0.7 fm^^ to 3.5 fm~^ is relatively small (1-2%), which is less than the uncertainty 
of the medium-effect parameter. In the case of the tin isotopes, the situation is similar to 
the case of ^^Ni. Thus, the effect from the coupled channel is masked by the uncertainties 
of our introduced medium effects in this analysis. Also, the coupled-channel effect might be 
partly included in our parametrization because we aimed to phenomenologically reproduce 
the experimental data. Thus, the effect of the coupled-channel is negligible and was not 
included explicitly in our analysis. 



IV. RESULTS AND DISCUSSION 



We have obtained the uncertainties of neutron density distributions considering statistics, 
systematic experimental errors, the uncertainties of the modification parameters, and our 
model, as shown in Fig. [91 Therefore the differences among the deduced shapes of the 
neutron density distributions of tin isotopes gradually changes. However the deduced RMS 
radii of the point proton and neutron density distributions clearly increase with the mass 
number, as shown in Fig. [TH This tendency is also shown in the results of several typical 
theoretical calculations, which are also plotted in the figure and have slopes consistent 
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FIG. 10: Coupled-channel effect of ^^Ni at 295 MeV. Solid lines show the coupled-channel calcu- 
lation between the ground state and the first 2^^ state. Dashed lines show the calculation without 
the coupled-channel effect. The global potential was used as the optical potential. Lower shows 
that ratios of cross sections and analyzing powers with to without performing the coupled-channel 
calculation. 

with our experimental results. Among the calculations, the nonrelativistic Skyrme Hartree 
Fock (SHF) calculations using SkM* parametrization are in good agreement with the RMS 
radii of both point proton and neutron density distributions. Neutron skin thicknesses Ar^p 
are shown in Fig. O The neutron skin thicknesses are about 0.11-0.19 fm for tin isotopes, 
which are not large what some RMF models suggest. The values of our Ar„p are reproduced 
by the SHF calculation using SkM* parameterization. On the other hand, the values are 
larger than SIII and smaller than NL3. Oyamatsu et al. and Chen et al. calculated a 
linear relation between the Ar„p and the symmetry term of the EOS 
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4l| . Therefore, 



our results favor medium values for the symmetry energy and its density dependence used 
in SkM* one, which are not so large as NL3 but larger than SIII. 
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FIG. 11: RMS radii of point proton and neutron of tin isotopes. Solid, dotted, and dashed lines 
show the results of theoretical calculations using typical mean-field models: RMF with NL3 
SHF including BCS force with SIII 



and with SkM* 



3lfi, 



371], respectively. 



V. SUMMARY 

In this work, we have extracted the neutron density distributions of tin isotopes. The ex- 
perimental data were analyzed in the framework of the RIA using a newly tuned relativistic 
Love-Franey interaction obtained for ^^Ni. The uncertainties of the neutron density distribu- 
tions were estimated in consideration of not only experimental but also model uncertainties. 
Using the tuned interactions in the nuclear medium, the neutron density distributions of 
the tin isotopes were deduced so as to reproduce the experimental data of the isotopes. 
We also deduced the RMS radii of the point proton and neutron density distributions. We 
compared our experimental results with those of several theoretical mean-field calculations. 
SHF calculations using SkM* parameterization as above were in good agreement with the 
RMS radii of both point proton and neutron density distributions. We observed a clear 
increase in neutron skin thickness with mass number, although the values obtained were not 
large what some RMF models suggest. 
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FIG. 12: Neutron skin thicknesses of tin isotopes obtained by various methods. Our results are 
indicated by sohd squares. Results from proton elastic scattering at 800 MeV [oj], giant dipole 



resonance 



38( 1 . spin dipole resonance 



39|], and antiprotonic x-ray data ||3$] are shown by open 
triangles, open diamonds, open crosses, open squares, respectively. The lines represent the models 
described in Fig. [TTl 
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